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Recent angle-resolved photoemission spectroscopy measurements of Co-doped LiFeAs report a 
large and robust superconducting gap on the F-centered hole band that lies 8 meV below the Fermi 
level. We show that, unlike a conventional superconductor described by BCS theory, a multiband 
system with strong interband Coulomb interactions can explain these observations. We model 
LiFeAs with a five-band model in which the shallow hole band is coupled with the other bands by 
only Coulomb interactions. Using Eliashberg theory, we hnd reasonable interaction parameters that 
reproduce the Tc and all five gaps of LiFeAs. The energy independence of the Coulomb interactions 
then ensures the robustness of the gap induced on the shallow band. Furthermore, due to the 
repulsive nature of the Coulomb interactions, the gap changes sign between the shallow band and 
the other hole pockets, corresponding to an unconventional s± gap symmetry. Unlike other families 
of iron-based superconductors, the gap symmetry of LiFeAs has not been ascertained experimentally. 

The experimental implications of this sign-changing state are discussed. 


I. INTRODUCTION 

Despite the differences between LiFeAs and other 
families of iron-based snperconductor^i^, the supercon¬ 
ductivity shares much in terms of family resemblanc^. 
Experiments have found that the anisotropic gaps in 
LiFeAs predominantly arise from antiferromagnetic spin 
ffuctuation^^HH These spin fluctuations originate from 
scattering between the electron pockets located at the 
M-point and the hole pockets at the F-point. In ad¬ 
dition, theoretical and first-principles calculations have 
confirmed the importance of spin fluctuations in LiFeAs, 
predicting an gap symmet^ similar to that of many 
iron-based superconductors^tlSI. 

However, the superconductivity mechanism at the in¬ 
nermost hole pocket at the F-point of LiFeAs remains 
a puzzle. Since the band barely crosses the Fermi level, 
the small size of the pocket makes studying it challeng¬ 
ing. Theoretical calculations have shown that spin fluc¬ 
tuations alone are insufficient to account for the large 
gap found on the tiny pockeiP. Furthermore, unlike 
the other pockets, the innermost hole pocket has an 
isotropic gajpp suggesting the presence of a different 
pairing mechanism. 

Recently, high-resolution angle-resolved photoemis¬ 
sion spectroscopy (ARPES) measurements have found 
that the gap on the innermost hole pocket is robust, 
even as the shallow band sinks 8 meV below the Fermi 
level upon electron dopin^SI. As shown in Figure 
the gap remains large, and the band’s spectral weight 
significantly changes between the normal and supercon¬ 
ducting states, even at energies well below the Fermi 
level. These observations suggest that the gap on the 
shallow band does not arise from low-energy excitations 
related to the structure of the Fermi surface. 

Prior to the ARPES measurement^^^ theoretical cal¬ 
culations have pointed to orbital-spin fiuctuationJi^ 
and the renormalization by high-energy excitationJi^ 
as means to obtain the large gap on the shallow band. 
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Figure 1. ARPES results of 3% Co-doped LiFeAs from Ref. 
El Left: the ARPES intensity plot in the normal state, 
showing that the shallow hole band at the P-point is 8 meV 
below the Fermi level. Middle: the same plot in the super¬ 
conducting state, showing that the large gap on the shallow 
band is robust against the Lifshitz transition. Right: a com¬ 
bined plot of the shallow band’s spectral weights in both 
the normal and superconducting states, showing significant 
changes between the two states even well below the Fermi 
level. 


However, it is unclear whether any of these proposed 
mechanisms will allow the large gap to be robust across 
the Lifshitz transition upon electron doping. 

In this paper, we propose Coulomb interactions as the 
superconductivity mechanism at the shallow hole band 
centered at the F-point in LiFeAs. We represent LiFeAs 
by a five-band model, in which the shallow band couples 
to the other bands by only Coulomb interactions. Us¬ 
ing Eliashberg theory, we find that interband Coulomb 
interactions induce a large superconducting gap on the 
shallow band. The energy independence of Coulomb in¬ 
teractions then ensures the robustness of the gap against 
changes in the Fermi level, in agreement with ARPES 
observationJlil Using reasonable interaction parame¬ 
ters, we find that our model can quantitatively repro¬ 
duce the experimental values of the and the gaps on 
all five bands. Relative to the case without Coulomb in¬ 
teractions, these interactions are found to enhance the 
Tc by a factor of 1.7, indicating the significant role they 
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play in LiFeAs. Finally, due to the repulsive nature of 
Coulomb interactions, our results predict an unconven¬ 
tional s± gap symmetry, in which the gap changes sign 
between the hole pockets at the F-point. Unlike other 
families of iron-based superconductors, the gap symme¬ 
try of LiFeAs has not been ascertained experimentally. 


II. METHODS 

Sinc e iron -based superconductors are moderately 
couplecpi^i3HI3 -^^0 employ Eliashberg theorjE^. Al¬ 
though our approach is phenomenological as in Ref. [9l 
BCS theory cannot be used, because it would be unable 
to yield the correct gap to Tc ratio. It is for this reason 
that we adopt the Eliashberg approach. 

In the Matsubara formalism, the multiband Eliash¬ 
berg gap equations ar^i^ 
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Here, i,j are the band indices, and ujn = (2n -|- 1) ttT 
is the nth fermionic Matsubara frequency. The function 
Zi{iujn) describes corrections to the electron self-energy, 
and Ai (icUn) is the energy-dependent superconducting 
gap of the Ah band. The potentials and de¬ 
fined by 
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describes the effects of band i on band j due to interac¬ 
tions with phonons and spin fluctuations, respectively. 
The Eliashberg functions a^Fij (w) can be experimen¬ 
tally determined from the inversion of tunneling data. 
In Eq. notice that the potential mediated by spin 
fluctuations appears with a negative sign. This is due 
to the spin-flip nature of magnon scattering. Finally, the 
pseudopotential describes the Coulomb interactions 
between bands i and j. It is usually given a cutoff at a 
large energy ojc for numerical convergence. Physically, 
the cutoff signifies the existence of an energy scale up to 
which the Coulomb interaction is instantaneous. Unlike 
boson-mediated interactions, the Coulomb interactions 
cannot be easily measured, and are usually tuned phe¬ 
nomenologically to fit experimental measurements. 

For a given set of interaction parameters, we self- 
consistently solve the Eliashberg equations for the 
superconducting gaps A^ (fw„). Then, using Fade 
approximants^Sl, we perform an analytic continuation to 
obtain the solutions in terms of real energies. The en¬ 
ergy gap Aq measured in ARPES experiments is then 
given by Aq = Re [A (w = Aq)]. 


III. FIVE-BAND MODEL 

We represent LiFeAs by a five-band model with a 
Fermi surface schematically shown in Figure The 
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Figure 2. A schematic of the Fermi surfaces of LiFeAs. In 
our model, the deep hole bands (1,2) are coupled with the 
deep electron bands (3,4) via interactions mediated by spin- 
fluctuations. In addition, the shallow hole band (5) is cou¬ 
pled with the inner deep hole band (1) via Coulomb interac¬ 
tions. 


deep bands ( 1 , 2 ) at the T-point are coupled to the deep 
bands (3,4) at the M-point by interband interactions 
mediated by spin-fluctuations, while the shallow band 
(5) at the F-point is coupled to the deep bands ( 1,2,3, 4) 
by Coulomb interactions. As in previous studies^J^ElIls]^ 
the spin fluctuations arise from interband scattering be¬ 
tween the hole and electron pockets. Phonon-mediated 
interactions involving only the deep bands are assumed 
to be negated by Coulomb interactions acting likewise. 
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and boson-mediated interactions involving the shallow 
band are assumed to be negligible due to the low density 
of states at the Fermi level. 

Previously, a four-band model of LiFeAs was stud¬ 
ied under similar assumption^. The authors omitted 
the shallow hole band claiming that its small density of 
states at the Fermi level precludes its contribution to 
superconductivity. They found that a large intraband 
phonon-mediated interaction on band 1 at the T-point 
is required to quantitatively reproduce the supercon¬ 
ducting gaps of LiFeAs. Here, we will show that this 
interaction is not necessary in our full five-band model. 

For the spin fluctuations in our model, we follow Ref. 
[5] and use Lorentzian Eliashberg functions (w) 

with peak energies 0^^ = 8 meV and half-widths Yij = 
4 meV. The coupling constants 
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satisfy A 31 /A 13 — 0.9019, A 41 /A 14 — 1.5010, A 32 /A 23 — 
1.0483, and A 42 /A 24 = 1.7447, in accordance with band 
structure calculation^. 

For the Coulomb interactions, since they are gener¬ 
ally stronger for smaller momentum transfers, we in¬ 
clude them only between the two innermost hole bands, 
that is only /ri 5 ,/r 5 i ^ 0. This can also be justified by 
the fact that only the two inner hole bands have similar 
orbital contentfA^l. Furthermore, intraband Coulomb in¬ 
teractions /Lt 55 can be omitted, as has been shown in a 
functional renormalization analyst^. While the values 
of can be calculated from first principles, doing so 
is difflcult, as they depend on the details of the band 
structure. The cutoff energy ujc of the Coulomb inter¬ 
actions is set to be lOOij, as is commonly done in the 
literatur^i^. 


IV. RESULTS 


Interaction 

parameters 

Energy gaps / 
meV 

n / K 

Ai3 

1.05 

Ai 

5.0 

Expt. 

18.0 

A 23 

0.73 

A 2 

2.6 

Model 

20.4 

Ai4 

0.41 

A 3 

-3.6 



A 24 

0.30 

A 4 

-2.9 



f^Sl 

0.38 

As 

-5.5 




Table I. The unique values of the interaction parameters 
Ai 3 , Ai 4 , A 23 , A 24 ,/ii 5 used to reproduce the five supercon¬ 
ducting gaps at low temperatures. The Tc resulting from 
this set of parameters is consistent with experimental mea¬ 
surements. 


with a, Tc ~ 20.4 that is consistent with experimental 
measurements. Such reproduction of the in addition 
to the gaps is often omitted in theoretical calculations. 
The larger errors at higher temperatures are expected, 
due to the reduced number of Matsubara points. Notice 
that the opposite signs between Ai, A 2 and A 3 , A 4 are 
due to the spin-flip nature of the interactions mediated 
by spin-fluctuations, while the opposite signs between 
Ai and A 5 are due to the repulsive nature of Coulomb 
interactions. 

To elucidate the effects of Coulomb interactions, we 
also performed low temperature calculations considering 
cases in which the Coulomb interactions /iy are scaled 
by a factor of 0 < a < 1. Figure shows that in the 
absence of Coulomb interactions at a = 0, the shallow 
hole band at the F-point is not gapped, A 5 = 0. As a 
increases, interband Coulomb interactions induce a gap 
on the shallow hole band, and increase the magnitude 
of the gaps on the other bands. For large enough a’s, 
the gap on the shallow band is the largest within the 
whole Brillouin zone. Then, because Coulomb interac¬ 
tions are energy independent in the regime of interest, 
varying only on the energy scale of the plasma frequency 
Wp ~ 1 e\“ the large gap on the shallow hole band is 
robust even as the band sinks below the Fermi level, in 
agreement with ARPES observationJSI. 

Next, Figure 1^ shows the effects of Coulomb interac¬ 
tions on Tc- Relative to the case without Coulomb in¬ 
teractions, these interactions are found to enhance the 
Tc by a factor of 1.7, indicating the significant role they 
play in LiFeAs. 


Although our model has six adjustable parameters 
Ai 3 , Ai 4 , A 23 , A 24 ,/ii 5 , and /rsi, we And that reproduc¬ 
ing the five gaps of LiFeAs requires pis > 0.2. We 
adopt the minimum value here to obtain a lower bound 
on the effects of Coulomb interactions. Now that the 
model is left with five adjustable parameters, they can 
be uniquely solved to reproduce the five superconduct¬ 
ing gaps at low temperatures. The results are shown in 
Tabled 

Using these parameters, we calculated the gaps at 
various temperatures, as shown in Figure The tem¬ 
perature dependence has the expected mean-field form 


V. DISCUSSION AND CONCLUSIONS 

While Coulomb interactions are often thought to op¬ 
pose superconductivity, our results illustrate the impor¬ 
tance of interband Coulomb interactions as a mecha¬ 
nism of Tc enhancement. Such enhancement is likely 
to occur in systems with bands that have small densi¬ 
ties of states at the Fermi level as Coulomb interactions 
dominate pairing for these bands. An example of such 
systems is single-layer FeSe, which has a hole band com¬ 
pletely below the Fermi leveP^. Pairing by interband 
Coulomb interactions may be important in the search 
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Figure 3. A plot of the energy gap against temperature for 
the five-band model of LiFeAs. The gaps agree with experi¬ 
mental measurements at low temperatures. The interaction 
parameters used are shown in Table [I] 



Figure 4. A plot of the energy gaps against the Coulomb 
interactions showing Coulomb interactions inducing a large 
gap on the shallow hole band. Due to the energy indepen¬ 
dence of Coulomb interactions, the induced gap is robust 
against changes in the Fermi energy. 
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Figure 5. A plot showing the effects of Coulomb interactions 
on Tc. In LiFeAs, Coulomb interactions enhance the Tc by 
a factor of 1.7. 


of new high-Tc superconductors. 

Our results also predict that the superconducting gap 
changes sign between the inner hole pockets at the F- 
point. The resulting sign-reversal gap symmetry is un¬ 
like the conventional s± gap symmetrj^^ believed to 
be present in iron-based superconductors. Similar un¬ 
conventional gap symmetries have also been proposed 
in Ref. UHl Unlike other iron-based superconductors 
in which the pairing symmetry has been extensively 
studiecP^, such is not the case for LiFeAs. Measurement 
of the gap symmetry can be performed using the SQUID 
junction proposed in Ref. [25l While there have been 
experimentgSSl measuring the gap symmetry of LiFeAs, 
the results are not conclusively, as the different hole 
bands were not individually resolved. 

ARPES experiment^y^ have also suggested that 
some iron-based superconductors are in the BCS-BEC 
(Bose-Einstein condensate) crossover regime, due to the 
large ratios of the superconducting gap to the chemi¬ 
cal potential observed. This observation can alterna¬ 
tively be understood in the context of our results. Since 
Coulomb interactions are energy independent, they can 
yield the observed superconducting gaps regardless of 
the size of the chemical potential, thereby providing a 
mechanism for the gapping of states lying below the 
Fermi energy. 

Recent experiments have provided further evidence 
that superconductivity in LiFeAs does not entirely arise 
from low-energy spin fluctuations. Nuclear magnetic 
resonance (NMR) measurements show that the Tc in Co¬ 
doped LiFeAs decreases even when the strength of spin- 
fluctuations increases upon dopin^^. Furthermore, tun¬ 
neling spectroscopy measurements show the existence of 
a temperature independent bosonic mode not directly 
related to spin fluctuation^^. These results strengthen 
our case that Coulomb interactions are important for 
superconductivity in LiFeAs. 

Before we conclude, we would like to highlight a com¬ 
mon misconception about multiband superconductivity 
found in the discussions of the ARPES results in Ref. 
m Unlike a one-band system, a gap driven by inter¬ 
band interactions in a multiband system depends on the 
density of states Nj of the other bands, and not on its 
own density of state W • For example, in a two-band sys¬ 
tem with only interband interactions, the ratio A 1 /A 2 
of the gaps is proportional to y/N 2 /Np^. Consequently, 
the shallow band in LiFeAs developing a large supercon¬ 
ducting gap does not contradict the principles of BCS 
theory, as was incorrectly implied in Ref. 1111 Neverthe¬ 
less, the paper’s main arguments remain sound. 

In conclusion, we proposed Coulomb interactions as 
the superconductivity mechanism at the shallow hole 
band centered at the F-point in LiFeAs. We represented 
LiFeAs by a five-band model, in which the shallow band 
couples to the other bands by only Coulomb interac¬ 
tions. Using Eliashberg theory, we found that interband 
Coulomb interactions can induce a large superconduct¬ 
ing gap on the hole band. The energy independence 
of Coulomb interactions then ensures the robustness of 
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the gap, in agreement with the ARPES observation^^]. 
Using reasonable interaction parameters, we found that 
our model can quantitatively reproduce the experimen¬ 
tal values of T^. and the gaps on all five bands. Relative 
to the case without Coulomb interactions, these inter¬ 
actions were found to enhance the by a factor of 
1.7, indicating the significant role they play in LiFeAs. 
Finally, due to the repulsive nature of Coulomb inter¬ 
actions, our results predict an unconventional s± gap 
symmetry, in which the gap changes sign between the 
hole pockets at the P-point. This study should help mo¬ 
tivate further experiments on the pairing symmetry of 


LiFeAs. 
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